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ABSTRACT 

Context. The observed light curves of most eclipsing binaries and stars with transiting planets can be well described and interpreted by 
current advanced physical models which also allow for the determination of many physical parameters of eclipsing systems. However, 
for several common practical tasks there is no need to know the detailed physics of a variable star, but only the shapes of their light 
curves or other phase curves. 

Aims. We present a set of phenomenological models for the light curves of eclipsing systems. 

Methods. We express the observed light curves of eclipsing binaries and stars, transited by their exoplanets orbiting in circular 
trajectories, by a sum of special, analytical, few-parameter functions that enable fitting their light curves with an accuracy of better 
than 1%. The proposed set of phenomenological models of eclipsing variable light curves were then tested on several real systems. For 
XY Bootis, we also compare in details the results obtained using our phenomenological modelling with those found using available 
physical models. 

Results. We demonstrate that the proposed phenomenological models of transiting exoplanet and eclipsing binary light curves applied 
to ground-based photometric observations yields results compatible with those obtained by the application of more complex physical 
models. 

Conclusions. The suggested phenomenological modelling appears useful to solve a number of common tasks in the field of eclipsing 
variable research. 

Key words, stars: variables - stars: eclipsing: - exoplanets 


1. Introduction 


There are two groups of astrophysical tasks standardly solved 
by the analysis of observations of eclipsing binaries (EBs) and 
stars with transiting extrasolar planets. We can derive the physics 
of the present state of a system and its components - namely 
the dimensions and geometry of the system, outer characteris¬ 
tics of both eclipsing bodies such as their radii, shapes, masses, 
temperatures, limb darkening, gravity darkening/brightening, 
gravitational lensing, albedos, spottiness, pulsation, and pa¬ 
rameters of possible streams and disks. Required information 
is extracted by means of various, more or less sophisticated 
physical models of double systems fe.g. [Wilson & Devinne ' 


physical models ot double systems (e.g. |Wilson & Devinne\ 

197 ll: iKallrath & MilonelfTo^ l Hadravall2004jBradstreel]l200-5l 




Prsa & Zwittej 20051: iPrsa et al. 1201 ll: Pribulld 20121) applied to 

numerous and precise data of all kinds obtained by contempo¬ 
rary observational methods and approaches developed for this 
purpose. 


Equally important are the description and classification of 
light curves (LCs), or studies of the evolution of the sys¬ 
tems, within the time scale of decades. These usually very 
tiny changes in light variation parameters (typically the period) 
may provide information on e.g. the rate of mass ex c hange 
between interacting co mponents (e.g. IZhu et all 1201 Ol 1201 2t 
iMikulasek et al.ll2012ah . the presence an d characteristics of pos¬ 
sible i nvisible bodies ( stars, planets, e.g. I Van Hamme & WilsonI 
llOO^IOian et^lIOOl in the system, or the degree of the mass 
concentration in stellar interiors (internal structure constants) if 
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we study apsidal motion in some eccentric EBs (e .g. lKopaJl978l: 
IClaret & Gimenene3ll993l : IZasche & Wollll2013h . etc. 

In period analyses of EBs we do not need all the informa¬ 
tion about the physics of the system, but only g ood templates of 
phase curves of all data we want to analyse (IMikulasek et al.l 
l2012al) . Such template phase curves (typically light curves 
in various photometric bands) can be obtained in principle 
by means of the sophisticated versions of physical EB mod- 
els applied to the complete set of observational data (s ee in 


IVan Hamme & Wilso 

nll2007l:IWilson & Va 

n Hamn 

iel20l4 

their best parts (e.g. 

Zasche et al.l 12014: 

Zaschel 

20L5h. 1 


is also the possibility of u sing the best observed LCs them¬ 
selves (iPribulla et al.l 1200 . We offer another alternative using 
relatively simple phenomenological (mathematical) modelling 
of observed data variations. 

Our aim is to establish a general model of light curves of 
eclipsing systems - ES (both eclipsing binaries and stars with 
transiting planets) that could fit the LCs with an accuracy of 1% 
of their amplitudes or better that can be applied to the major¬ 
ity of obse rved eclipsing systems (the concept of the model is 
outlined in IMikulasek et al]l2015l) . Such a model could be used 
for most LC description tasks and detailed period analysis, for 
which other sources of phase information could be also used, 
especially individual eclipse timings and radial-velocity (RV) 
curves. 

The remainder of this paper is organised as follows: in Sect.|2] 
we describe the general properties of phase variation of pe¬ 
riodically variable objects, Sect.|3] specifies general properties 
of eclipsing system LC models, and Sect.|4] presents the phe¬ 
nomenological model of eclipsing system light curves. Sect.|5] 
is devoted to model parameters and searching for them dur- 
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ing analysis, in Sect.|6] we compare the results of the mod¬ 
elling XY Boo ob servations processed by mean s phenomenolog¬ 
ical and physical (IWilson & Van HammdlZOl^ methods, and in 
Sect.|2]we summarise and discuss our results of the phenomeno¬ 
logical modelling of the eclipsing systems. 


2. Periodically variable objects 


Eclipsing binaries and stars with transiting exoplanets periodi¬ 
cally change their light mainly due to regular mutual eclipses 
of components (transits and occultations) and proximity effects 
(tidally induced ellipticity of components and reflection). The 
period of these ES variations corresponds to the observed orbital 
velocity of the system; the shape of light curves, dominated by 
eclipses and proximity effects, is hence more or less constant 
(for details see Sect.O. That is why ESs are ranked into a group 
of periodically variable stars, or objects in general. 

Most of the variations of periodic variables are more or less 
cyclic with an observed instantaneous period Pit) which is usu¬ 
ally strictly constant or slightly variable. The period itself and its 
development over time can not be observed directly, but both can 
be derived through the analysis of time series of light changes or 
extremum timings. Eor that purposes it was useful to introduce a 
monotonically rising phase function ??(f), as a sum of the epoch 
E{t) and t he phase in the commo n usage <f> and its inversion func¬ 
tion t{ff) (iMikulasek et alJl2008l) . 

Eunctions d-{t) are tied with the instantaneous observed 
period P(t), P(§) by the following differential eq uations (see 
iKalimeris et alJfT^ l^ulasek et al.ll2008ll2012bl) : 


d)?(f) 

df 


1 

m’ 


d{Mo) = 0; 


dt(-&) 

d§ 


= p(py, m = Mo, 


( 1 ) 


where Mq is the origin of counting of epochs (for ESs the time 
of the basic primary minimum). 

Using tiP) we can predict the zeroth phase time 0{E) for the 
epoch E (primary minimum timing according to the ephemeris 
Eq. O), 0{E) = 0(P = E). 


2.1. Period models. Phase and time shifts. 


The linear period model supposes the period P(t) of the variable 
object to be constant. The cotTesponding linear phase function 
i?i(f, Mo, Po) and its inversion ti{§, Mo, Pq) are then given by: 

§iit) = it- Mo)/Po, ti{§) ^ Mo+ Po §. (2) 


A possible tiny modulation AP(f) of the basic period Pq causes 
a detectable phase function shift Ad(t) in light curves and shifts 
At(P) - A0{E) in EC extrema timing. 


Pit) = Po + AP(f); Pit) ^Py+ Atpit); tiP) = ti+ AtiP). (3) 


Combining definitions in Eq.j^with the fact that functions P(f) 
and tiP) are mutually inverse, we obtain the following relations 


At = -Po AP; 

APit)± f 
Jmo 


APiP) = 
AP(t) 


dAf(i?) 


d)? 
(AP(t))2 


p2 


p3 


d[f(P)-ti(7?)] 

dP 

AP(t) 


d.. f 

Jmo 


pi 


dr. 


(4) 

(5) 


The last approximation in the Eq. (|5]l is valid for all known EBs, 
including SV Cen with a record-breaki ng decrease of its orbi tal 
period by P/Po = -2.36(5) x 10“^ yr * (IDrechsel et alJll98^ . 

The time development of phase function shifts APf) can be 
derived by the analysis of light curves, whilst A6>(P), can be 
revealed using standard 0-C diagrams constructed by means of 
extrema timings. 


2.2. Basic period models. Finding of 0-C shifts 


If the time development of the period Pit) is continuous and 
smooth we can express it in the form of the Taylor polynomial 
with the centre at f = Mq, Pi - it - Mo)/Po 


Pit) = Po + PoPo + EgEo J + ... + ^0^1 ■ • • (6) 

Using Eqs. (I2]l, Q, (|4]l, and simplified version of Eq. (|5]l we 
obtain: 


^k-\ dTo 




§ ^ Pq 2 , PqPq 31 ... Pq 


t ^ Mq+PqP + PqPq^+pIPq^ + . . 

ok dTo 


0 ^ Mq+PqE + PqPq^+P^Pq^ + . . ■+f^o“d(^ (I+T)T 


(7) 

( 8 ) 


Similarly we can establish other arbitrarily complex period mod¬ 
els of Pit) determined by a set of free parameters including 
also cyclic period modulation of the light time effect (LiTE) 
caused by anothe r body in the system (IMikulasek et al.lf201 Ibl: 
iLiska et d]l2015l) or the apsidal motion. 

The real shape of the phase curve can also be approximated 
using so-called 0-C time shifts of the observed phase curves ver¬ 
sus the predicted light curve (EC) derived by the period model 
with fixed parameters (typically Pq, Mq) expressed in time units 
(usually in days). Let us divide the whole time interval covered 
by observations into noc appropriate time intervals (typically 
nights or seasons). The phase function Pit, r) during a certain 
r-th time interval with the jO-O r time shift is then given by the 
formula (IMikulasek et al.ll201 lah 


iO-Qs 

Pit, r, P,, {0-C}) = PM - 2] rjrs (9) 

.?=! 


where Pdti) is a predicted phase function calculated by an ap¬ 
propriate period model at the time f,, [O-C] is a set of all noc 
values of the found iO-C)r time shifts versus this model. The 
symbol rjrs represents a discrete function (a table or a matrix), 
which for each individual i-th observation from our data assigns 
either 1, if its order number of the interval r, = s, or 0, if r, s. 

The dependence of O-C values on the epoch serves as a 
common O-C diagram, the basic tool for the period analysis. 
Using the found individual O-C values and their uncertain¬ 
ties we can calculate a set of so-called virtual minima tim¬ 
ings (IMikulaseket al.l 1201 lafl or Sect. 16.4.21 in this paper). 
Virtual minima timings can be combined with others derived 


e.g. by other techniques (e. 

g. Brat et al. 20121 Mandel & Agoll 

2002, 

Kwee & van Woerden 

119561 Mikulasek et al .1120061 12014; 

Zasche et al.ll2014l: Zaschell2015l). 


3. General properties of an eclipsing system light 
curve model 

3.1. Instrumental term of an observed light curve 

The observed light curve (or its segment) of a chosen eclipsing 
system in a particular colour of an effective wavelength Teff is 
defined by a time series {f^,, y„} obtained during an observational 
interval r (observing night, part of it or season). The r-th subset 
of observational data can be generally modelled by the function 
Yrit, deff), expressed in magnitude^ 

YAt, deff) = YoAt, deff) + FiP, deff); (10) 

* The presented models may also be applied for expressing light in¬ 
tensity variations, but the treatment of data in the magnitude domain is 
more straightforward and the accuracies of results are almost the same. 
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where ForCf, ^etr) is an additive term removing instrumental 
trends and concurring the observed magnitudes or magnitude 
differences {y,}, whilst F(§, /left) is an intrinsic light curve model 
function of the phase function ■& (see Sect. l2.1l i and the effec¬ 
tive wavelength deff, free of instrumental and observational shifts 
and trends. The function F(j?, Teff) is normalized so that its mean 
value without eclipses is equal to zero. The function Yor(t, Tefr) 
can be approximated by a linear combination of gr dimension¬ 
less functions of time t, Ej (t), with magnitude-like coefficients 

gr 

Tor(f»Teff) — ^ JflQi-j {t), (11) 

1=0 

where Ej (f) can be e.g. normalised polynomials Eyj (f) = [(f - 
ir)lo'{tr)V {cr{tr) is the weighted variance of observational times 
in the segment r) or Legendre polynomials or special quasi or¬ 
thogonal fun ctions combining polynom ials with harmonic func¬ 
tions (see in iMikulasek & Gril l2005h . The set of coefficients 
{m^jrj (Teff)} is found together with parameters describing the 
model function F{§, Teff). 



Fig. 1. Examples of elementary model LC W = for C = 0, 

see Eq. GS allowed by our phenomenological EB models. Parameter 
d describes the duration of the eclipse, T is a parameter expressing the 
kurtosis of the LC. 


3.2. Bases and limitations of the phenomenological model of 
eclipsing system light curves 


Light curves of ESs are nearly periodic functions and it 
would be nat ural to express them in the form of Fou r ier se¬ 
ries ( s ee e.g. Rucin^ [l^T^ lKullrath_& Miltm j 1 19981 : ISelamI 

ll 


l2004t iNedOTOscrk et al.l |2012t Andronov 20121) . This concept 
proves its worth in many types of extrinsic periodically vari¬ 
able stars, especially in the case of rotating variables with pho¬ 
tomet ric spots on their surfaces dNorthl 11984 IMikulasek et al.l 
l2007l) or non-eclipsing (e.g. purely elliptical) double stars, where 
we manage with harmo nic polynomials of a low degree (e.g. 
iKallrath & Milonelll998h . 


However, it is generally known that the presence of eclipses 
in LCs asks for the use of harmonic polynomials of a rather high 
degree if a proper fit of the observed LC is required. Their us¬ 
age is badly influenced by departures from the ideal equidistant 
distribution of observations according to the orbital phase. Even 
small phase gaps are then filled with unreal LC artifacts. 

It seems that it is better to use properly selected phenomeno¬ 
logical models of LCs of eclipsing systems described by a 
few parameters. Several more or less successful at tempts on 
how to model these LCs have been proposed e.g. by iTsesevichl 
(ll97lh:lKholonovl (119811). see also the re views and references in 
lAndronovl(l2Q12l) : IChrastina et al.l (l2013h . Table|2] and Fig.[T]|3 


Theory is able to explain the observed periodic light vari¬ 
ations of an ES in general and in detail as the result of al¬ 
ternating mutual eclipses of components of the system, non¬ 
isotropic radiation of orbiting components, caused by their 
close proximity, and by uneve nly distributed photometric spots 
on rotating compon ents (e.g. iBradstreel 120051: iHilditchIbOOll: 
iPribulla et al.l l2Ql^ . The periodicity of light changes caused 
by mutual eclipses and proximity effects, as well as variations 
connected with the rotation of synchronously rotating spotted 
components is dictated by their observed orbital motion. These 
changes are periodic with an instantaneous periocQ P(t), the 
shapes of light curves remain more or less constant for decades. 

For purely geometrical reasons the prevailing majority of 
EBs ranks among relatively close, and hence tidally interacting 


^ Observed orbital period may change due to possible transfer of mat¬ 
ter between components or light time effect. 


systems, where the processes of synchronizing of the compo- 
nents’rotation and orbit circularization are strong and effective 
(IZahnlll9^ iGoldman & Mazehlll99lL and references therein). 
That is why the rotations of components are usually synchronous 
and mor e than 80% of their orbits are pretty circular (see 
CALEB. iBradstreei 1200-51) . In the further introductory text we 
will concentrate mainly on systems with more or less circular 
orbit^ and constant light curvefl 

The intrinsic one-colour light curve function (expressed 
in magnitudes) of an eclipsing system is a periodic function 
F(j 7, Teff) that can be approximated as the sum of th ree more 
or less independent terms (see also in I Andronovl20 1 2|) : 


Fi§, deft) = F,(&, deft) + Fp(P, deft) + F,(P, deft), (12) 


where Fe(&, deft) describes a light curve of eclipses (Fg = 0 out¬ 
side of eclipses), whilst deff) and Fe(f?, d eft) express con¬ 
tributions of proximity and lO’ConnelH (Il979l) effects without 
eclipses (Fp - Fg = 0). The mathematical models are formu¬ 
lated and discussed in Sect. l4.L2l 


^ Applying so-called phase rectification (Mikulasek et al. in prepara¬ 
tion) we are also able to solve eccentric systems. The technique sym¬ 
metrizes light curves of eclipsing binaries with components moving un¬ 
evenly on their eccentric orbits by the rectification of their phase func¬ 
tions. The method of phase rectification enables also the effective anal¬ 
ysis of apsidal EB motion. 

This assumption is fulfilled only partially. We can mention cycli¬ 
cal variations of instantaneous LCs with other than orbital period as the 
gradual change of the ES geometry due to asymptotic motion of dou¬ 
ble systems orbiting on eccentric orbits, possible asynchronous rotation 
of spotted components in wide systems, and possible pulsations of the 
components. Eclipsing binary light curves may also vary erratically be¬ 
cause of chromospheric activity (see e.g. Sect l6.4.^ . time-dependent 
spottiness of the components, or changes in streams or disks around the 
stars. Neglecting of above mentioned effects introduce as a rule some 
extra noise in the period analyses and deteriorate the accuracy of the 
determination model parameters. 
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phase 


Fig. 2. The fit of the V curve of the overcontact spotted EB EKCom 
(the secondary minimum is a transit). The LC of the close binary is 
affected by O’Connell and proximity effects. The synthetic model curve 
is depicted by a solid line, the phenomenological LC fit is depicted by 
the dashed line. Eor more information see Sect. l4.3l and Table[T] 



Fig. 3. The simultaneous fit to 15 transits of extrasolar planet TrES- 
3b corrected for trends (see Eq. dlOH . The fit parameters according to 
Eq.dc A = 0.02725(19) mag, D = 0.0117(7), C = -0.17(7), and 
r = 1.58(12), the parameter K was not introduced. The black line is 
the fit, circles are the normal points, grey circles indicate individual 
measurements with the area proportional to their weights. 


4. Phenomenological model of eclipsing system 
light curves 

4.1. Model of one-colour light curves 

4.1.1. Eclipses 

The essential feature of all ES light curves are two nearly 
symmetrical depressions caused by mutual eclipses of syn¬ 
chronously rotating stellar or planetary components. The pro¬ 
files of both minima are complex functions determined primar¬ 
ily by the geometry of the system and the relative brightness of 
components in a given spectral region centered at the effective 
wavelength deff. The contribution of eclipses Fe()?, deff) to an ES 
light curve can be approximated by a sum of two special peri¬ 
odic functions of phase function §. In the case of circular orbits 
eclipses are exactly symmetrical around their centres at phases 
1^01 and (fo 2 . If we put the origin of the phase function Mq at the 
time of the primary minimum, then (foi - 0, (po 2 = 0.5. 

The model function was selected so that it describes as aptly 
as possible those parts of LCs that are in the vicinity of their 
inflex points, where their slopes are maximal. The functions are 
parameterised by their widths Di, D 2 , eclipse LC kurtosis coef¬ 
ficients El, E 2 , dimensionless correcting factors Ci, C 2 , and cen¬ 
tral depths Ai(defF), A 2 (deff): 


WdefF)=£Ajl+a- 


k=l 




1 -•! 1 -exp 


1-cosh 


Dk 


<Pk ^ §-0.5 (k-l)- round [& - 0.5 (k - 1)], 



where the summation is over the number of eclipses during one 
cycle, n^: Ue = 2 or «e = 1 (the common situation for exoplanet 
transits). Each eclipse in a given colour is thus described by only 
four parameters - its depth A, width D, kurtosis F, and the cor¬ 
recting parameter C. 

In the case of eclipsing binaries with two minima in a cy¬ 
cle (Ue = 2) we need eight parameters, but sometimes the num¬ 
ber of needed parameters can be smaller. Inspecting the pa¬ 
rameters D, r, and C for both eclipses of many EBs we have 
concluded that they are as a rule nearly the same: especially 
Di = D 2 , Fi = F 2 , and Ci = C 2 . So we usually need only five 


monochromatic parameters (Ai, A 2 , Z), F, C). The parameter C is 
mostly comparable with its uncertainty, so we can neglect it en¬ 
tirely. Then we need just four parameters! On the other hand, in 
EBs with totalities we see that the bottoms of their occultations 
are flat whilst transits are convex. It can be described by intro¬ 
ducing of different parameters Ci, C 2 (see the case of EK Com 
in Tabled Eig.|2li. 

The LCs of the exoplanet transits (wg = 1) need only 
four parameters {A,D,F,C - see Eig.|3]l, in cases of very pre- 
cise measurements we add another dimensionless parameter K 
(iMikulasek et al.ir2015l) : 


Ee = A 



1 - exp 



. (14) 


Testing several dozens of LCs of various types of eclipsing sys¬ 
tems we found that the standard deviation of the fit is typically 
well bellow one per cent. The only minor inconvenience is the 
existence of a spike (a jump in derivatives) in the mid-eclipses 
for LCs with E < 1 (see Eig.[T]l. 


4.1.2. Proximity effects. O’Connell effect 

Light variations of EBs caused by eclipses are usually modi¬ 
fied by asphericity of the components, effects of gravity dark¬ 
ening/brightening and mutual irradiation. All these proximity ef¬ 
fects are the manifestation of the interaction between the compo¬ 
nents acting, namely, in close systems. Contrary to eclipses the 
proximity effects modify LCs permanently, in each phase. 

Light curves of some eclipsing binaries are influenced by 
the O’Connell effect that results in the asymmetry of some 
LCs of close EBs, manifesting as the difference in light max¬ 
ima between eclipses. The standard explanation for this is 
the presence of one or more cool or hot spots on the sur¬ 
face of one of the synchronously rotating components or by 
asy mmetrically distributed c i rcumstellar material in the system 
(e.g. IWilsey & Beakyl 120091 : iPribulla et al.l 1201 21 and citations 
therein). T he amount and the si gn of the O’Connell effect yary 
with time dBeaky & Koiull2()T^ . The effect is also wayelength 
dependent - in blue it is usually stronger, but it is not a rule 
(IPribulla et al.ll2012h . 
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Table 1. Parameters describing LCs of several eclipsing binaries. 


Name 

ARAur 

EK Com 

AVDel 

477 Lyr 


9 

10 

13 

8 

filter 

BjV 

V 

BVRl 

V 

n 

686/687 

247 

545 

294 

D 

0.0168(3) 

0.052(3) 

0.0396(7) 

0.0219(5) 

r 

0.619(10) 

1.7(4) 

0.82(5) 

0.79(4) 

an 

0.691/0.681 

0.331(12) 

1.153(11) 

1.557(14) 

an 

- 

- 

1.17(4) 

- 

an 

0.551/0.564 

0.289(9) 

0.123(3) 

0.078(5) 

an 

- 

- 

-0.188(13) 

- 

an 

-0.0022(5) 

0.023(5) 

- 

0.321(2) 

an 

-0.022(4) 

- 

- 

_ 

^41 

0.0028(9) 

0.150(5) 

0.050(2) 

-0.035(2) 

(242 

- 

- 

-0.044(1) 

- 

(251 

- 

- 

- 

-0.005(2) 

(261 

- 

0.032(2) 

- 

- 

c 

-0.25(2) 

- 

-0.24(7) 

-0.342(20) 

Cl 

- 

0.10(17) 

- 

- 

Cj 

- 

-0.14(9) 

- 

- 

p 

1.05 

0.73 

1.05 

0.31 

type 

OC 

SD 

D 

D 

P 

4“; 13169 

0;26669 

3‘;8534 

0147173 

fi 

0.098 

0.524 

0.198 

0.077 

ri 

0.100 

0.318 

0.343 

0.219 

Teffl 

11 100K 

5000 K 

6000 K 

60 000 K 

Tefn 

10 600 K 

5300 K 

4275 K 

6500 K 

i 

88.5° 

88.5° 

81.3° 

65.8° 

Fig. 

13 

m 

m 

i 


ga is the number of the parameters used for the description of a LC, p = 
dphen/«5phys. where (5phen ^nd (5phys are the uncertainties in the determination of 
the zero phase time according to phenomenological and physical models. The 
used‘hi-fi’ LC model: A) = 2^1 [1-t C* ((/:*/D)2] 11 - {1- 

exp [1 - cosh(^ot/D)] 1^1 + 2^=3 cos [2'K(k-2)-ff\+aija sin(2ff i?). 

EB types: OC - overcontact, C - contact, SD - semidetached, D - detached, P is 
the period, ri,r 2 are the relative radii of the components, Teffi, 7’eff2, their ef¬ 
fective temperature, i is the orbit inclination. The parameters were taken from 
authors cited in Sect. 1431 

The contribution of proximity effects should be 

an even function symmetric with the phases 0.0 and 0.5 
and consequently they can be satisfactorily well expressed 
as a linear combination of elementary cosine functions 
cos(27r j?), cos(4;7r)?), cos(6;7r)?),... . The even terms are the 
consequence of the ellipticity of tidally interacting components, 
whilst the odd terms result from the differences between the near 
and far sides of components. As a rule we c an limit ourselves 
only t o the first two or three term s in the dRussell & Merrilll 
IT9^ iKallrath & Milonel Il998h . The O’Connell effect con- 
tribution Friff) can be well modelled by a si mple sinusoid 
(iDavidge & Milonel 1984tlWilsev & Beakvll200^ : 

«p+«e «c+Wp+«e 

F-p - ^ Aj: cos [2;7r(A:-Ue))?'], ^ A*: sin(27rj9'), (15) 

/:=«e+l ^=«p+«e+l 

where tip is the number of terms in Fp{§): tip - 0,1,2,3,..., 
Wc = 0, if the O’Connell asymmetry is not present, else Wc = 1- 
The V light curve of the close eclipsing binary EK Com (see 
Fig.|2]) with apparent O’Connell effect is described by nine pa¬ 
rameters: Ai, 2 , 3 , 4 , 5 ,£>, Ci_ 2 ,C, tie - 2, tip = 2, tic = 1. The un¬ 
common V light curve of EB V477 Lyr (see Eig.|4]l is determined 
by eight parameters: Ay 2 , 3 , 4 , 5 , D, C, F, ti^ - 2,tip - 3, ti^ - 0. 

^ If p > ^ then YjI^p hk = 0. 



Fig. 4. V light curve of a famous EB consisting of a very hot nucleus 
of a planetary nebula and solar type star. For the fit of the LC strongly 
affected by proximity effects we need only 8 parameters phenomeno¬ 
logical. For details see Sect. l4.3l and Table[T] 



phase 


Fig. 5. The fit of the primary minimum of the detached EB AR Aur in 
V and B, represented by normal points (circles), by physical (grey solid 
lines) and phenomenological (dashed lines) LCs (see Table[T](. 


4.2. The model of multicolour light curves 

The parameters of the above defined model functions, especially 
the amplitudes A^, and parameters Cip, F>i 2 and Fi ^2 are gener¬ 
ally functions of the wavelength A. 

In principle, we can use the one-colour models formulated 
in Sect. l4.ll separately, assuming that all of the parameters are 
wavelength dependent. Fortunately, it follows from our experi¬ 
ence with modelling of LCs of hundreds of real systems and 
their physical models that it is not necessary to take response 
curves of different photometric passbands into account; we man¬ 
age with their effective wavelengths only. It enables the associ¬ 
ation of photometric colours with different transparency widths, 
equal or close effective wavelengths (typically we are allowed to 
combine measurements done in V and y). 

On top of that the dependencies of model parameters on the 
effective wavelength Teff are typically smooth, mostly mono¬ 
tonic, so we can approximate them by low-order polynomi¬ 
als of the dimensionless parameter A; A = do/deff - 1, where 
do is an arbitrarily selected central wavelength of the data set 
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Fig. 6. Ten- colour transit LCs for HD 209458b taken from 
iKnutson et alJ ll2007h . Effective wavelengths of individual colours 
are in nm. For the desc ription of all 10 light curves we need only 9 
parameters (courtessy of lMikulasek et al.ll2015h 



Fig. 7. The dependence of the amplitude A of transit of HD 209458b 
(see Fig.[6l( on A = 550//leff - 1 proves the validity of the approximation 
Eq.GU for that case. 


(iMikulasek et al.ll2015l see also Fig.|7]i: 

= Z%\ akjAj-\ Cl = Y.% (16) 

Di^Y,%\dijAi-\ r,^Yf[liyijAi-\ 

where gAk, ga, goi, gn, I = 1,2 or / = 1, are the numbers 
of degrees of freedom of the corresponding parameters of the 
model. The standard set of the one-colour LC model parame¬ 
ters of EBs (see Sect. l4.11 l: {A^, Ci,Di,ri] can be considered as 
the special case of the multicolour decomposition Eq. (fThl) for 
do = deff, A -0, j - 1: {aki,cn,dn,yn}- 



AV Del 

- 0.2 0 0.2 0.4 0.6 

orbital phase 

Fig. 8. The fit of the BVRJc LCs of a semidetached short periodic 
‘cool Algol’ AV Del. Dark dots are observations, grey lines is the 13- 
parameter phenomenological fit. For details see Table[T]and Sect. l4.3l 


The set of relations in Eq. (fThl) enables the calculations of all 
parameters {Ak{A), Ci(A), A(A), //(A)) needed for calculation of 
the model of a LC in any photometric band characterized by the 
parameter A. 

Fig -0 shows the fit of the extreme 10 -colour photome try 
(350-980 nm) of an exoplanet transit dKnutson et al.l 1200^ - 
here we need 9 parameters, namely A - ai + a 2 A, D - di+ d 2 A, 
C - C] + cpA -H c^A^, r - y\ + 72 A, where A = 550/deff - 1 
dMikulasek et alfcOlSl) . 

The fit of BWch proper LCs of AV Del (see Eig.[^ needs 
only nine parameters, namely D\ - D 2 - D, Fx - r 2 - F, Ci = 
C 2 = C, r^ii, kii 2 i ^ 21 , ^ 22 , ^41, and ^ 242 , ne=2, iip — 2,Uc = 0, 
A 3 s 0. See also Table[T] 

4.3. Brief description of seiected eciipsing systems 

TrES-3b is an extrasolar planet orbiting the star GSC 03089- 
00929 with a period of 31 hours. It belongs to the hot Jupiters 
which are undergoing orbital decay due to tidal effects. Eor the 
LC inspection (Eig.O R photometry of 15 transit s containing 
2820 individual measurements in total (courtesy IVanko et’aP 
l2013h was used. The parameters of the LC fit of TrES-3b are 
in the legend of Eig.[3 

HD 209458 was the fir s t star found to have a transiting planet 
dCharbonneau et aP I 2 OOOI : iHenrv et~^ I 2 OOOI 1 and remains the 
second brightest star known t o have a transiting planetary com¬ 
panion. l^utsDn^eT^ d2007h obtained 1066 spectra over four 
distinct transits with the STIS spectrometer on HST allowing 
to synthesize LCs in 10 spectrophoto metric bandpasses in 29 0- 
1030 nm (see also Sect. l4t21 Eig.0and lMikulasek et al.ir2015h . 

Tab le[T] contains the parameters of the phenomenological fit 
of LCs and some other information on the following selected 
eclipsing binaries: 

0 AR A urigae, a prototype of a detached EB, dO’Connelll 
Il979h . 

o EK Co mae, an overcontact, spotted EB with a short orbital 
period, dSamec. Gray & CarriganllTo^ . 

o AV Delphinis, a ‘cool Algol’ consisting of a E type primary 
on the main sequen ce and a K subgiant filling its Roche lobe 
dMaderetal.ll20n.5h . 
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coefficients (typically the fits of LCs or minima times) should 
be computed by the general law of uncertainty propagation as- 
suming also correlations amo n g individual coefficients (s ee e.g. 
iBevington & RobinsonI l2003t Mikulasek & Zeid"^ l2013h . It is 
advisable to orthogonalize the models at least in the ephemeris 
parameters, in accorda nce with what we did in (lMikulasekl2007l: 
Mikulasek et al.ll2()()^ . 

5.2. The selection of an optimal model of light curves 


Table 2. The list of alternatives for several phenomenological model 
functions of eclipses. 


Fig. 9. The difference in mag between the physical model primary min¬ 
imum of ARAur and various alternatives of phenomenological basic 
function listed in Tabled 


o V477 Lyrae, an unusual, detached EB consisting of a very 
hot and luminous nucleus of the planetary nebula as a 
primary component a nd a solar type star as secondary 
(iPollacco & Belllll994 . 

Table[T] shows that the fit of LCs of the above mentioned stars 
by our ‘hi-fi’ models, quantified by the ratio p is nearly the same 
or bett er than in the cas e of the fit of BM3 (CALEB) physical 
model (lBradstreeljl2005h . 

5. Phenomenological model solution 

5. 1. Finding of model parameters and their uncertainties 

The procedure for finding model parameters is based on the si¬ 
multaneous mathematical processing of all relevant photometric 
data consisting of individual photometric observations, includ¬ 
ing barycentric Julian date of the /-th measurement f,, the mea¬ 
sured magnitude or magnitude difference corrected for possible 
trend(s) during nights or seasons y,, and the estimate of its un¬ 
certainty CT/. Lurthermore, we should know the effective wave¬ 
lengths of the photometric filter used, or A, , and submission 
of an individual observation to one of the observational subsets 
r, (see Sect l3.1l l. 

Lor simplicity we shall assume that the shapes of LCs are 
constant and the variability of an object is described by the 
unique model function Eq. (fTOl i. consisting of the instrumental 
term Yroft, Teff) (Eq. (fTTb l and the intrinsic phenomenological 
ES model EC function F(§, /left), specified in Sect.|4] The phase 
function §(t, Pq, Mq, Pq, or {0-C },...) is a function of time and 
some free parameters of a variety of period models offered in 
Sect.|2] The result of the solution - the full set of g free pa¬ 
rameters of the complete model including the estimate of the 
parameter uncertainties was evaluated simultaneously using the 
non-linear least square method by minimising the quantity by 
the standard technique (using tried Newton-Raphso n method of 
non-li n ear equat i on so l ution) well described i n e.g. IP ress et al. 
(20()2|): lHavashil (120001) : iHartkopf et°^(ll989l) : Mikulasek et al. 
( 201 Ibl) . With a good initial estimate of the parameter vectors 
the iterations converge fairly quickly. 

All estimates of uncertainties of model parameters were 
computed using formulae taking into account that our models 
fit phase curves of EBs with uneven accuracy. Since the mod¬ 
els are not orthogonal, uncertainties of the functions of model 



Fciv) 

s parameters 

-1-7TTTTT 
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5.5 

A,D,r 

6.9 

A,D,r 

11 

A,D 

16 

A,D 

18 

A,D 


The fidelity of individual alternatives are illustrated by the Fig.|2]and quantified 
by the scatter s of the fit of the ‘real’ light curve of primary minimum of AR Aur, 
simulated by its physica l model. The model function number 4 was suggested by 
I Andronov! 120101 the function ‘Real’ extracts a real part of the argument. 


The models of light curves should be tailored to the studied 
object/objects, available data, and the purpose of fitting the data. 
Especially, the number of used free parameters should be as few 
as possible, however without any serious influence on the accu¬ 
racy and reliability of the results. It follows from our experience 
that it pays off to adhere to the following general principles; 

1. It is always advantageous to process all available data si¬ 
multaneously and not divided into parts. It is also valid for 
the determination of mid-eclipse times of individual eclipses 
where we should use the Eq. (|9]l. 

2. We have to pay attention to the right weighting of entered 
data as it is required by the used regression. If our origi¬ 
nal data do not content individual uncertainties, we have to 
estimate them iteratively from the scatter of residuals {Ay,} 
for appropriately created data subsets. 

3. Eixing of EC model parameters, if they are known from pre¬ 
vious analyses, is also advised. In addition, neglecting of all 
model terms whose amplitudes are less or comparable with 
their uncertainties, is recommended. 

4. The number of free parameters can also be reduced by us¬ 
ing simpler model functions of eclipses listed in Table|2] 
Suggested model functions are compared in Eig.|9] 

5. The complexity of the selected phenomenological model 
should correspond to the purpose of the modelling. Several 
tasks (e.g. the basic classification of LCs) allow for the use 
of simple, unified models with basic parametric outfit. 

Unfortunately, the effort of using the optimum ‘hi-fi’ phe¬ 
nomenological models considerably encumbers automation of 
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the computational process. The diversity of real light curves of 
particular ESs requires that they be solved individually. 


6. Period analysis of XY Bootis in 1955-2009 


The above described phenomenological modelling of periodic 
variable stars has been developed first and foremost for the sake 
of the simultaneous period analysis of all available data contain¬ 
ing phase information. Such data are in the case of ESs stan¬ 
dardly times of light minima, LCs or their segments and radial 
velocity curves. The majority of studies of ES period changes 
were based on the analysis of 0-C diagrams constructed us¬ 
ing the times of light minima, where these timings were de¬ 
termined individually directly by observers not taking into ac¬ 
count LCs of the star observed before. Several other studies were 
based only on the analysis of LCs obtained during several years. 
Measurements of radial velocities were used, namely, for the so¬ 
lution of the geometry of the system. 

The simultaneous analysis of data of various kind is pos¬ 
sible using an extended version of the contemporary so¬ 
phisticated cod es for the phys i cal solution of eclipsing sys¬ 
tems (see e.g. iH adraval | 2004 IVan Hamme & Wilson 


I Wilson & Van Hammel 2014t) . The concept of EB period anal¬ 
yses without using 0-C diagrams a nd the physical s olution 
of the system was outlined (i.a.) in iMikulasek et al.l (l201 lai 
l2012ai 2013al). and studies of period changes of individual EBs 
dZhu et al.ll20^ I2012I: iMikuMsek et ani2013bh . The results of 
the preliminary report on the ephemeris of AR Aur - a star with 
the light time effect (IMikulasek et al. 201 Ibh - were r eferre d 
to and extensively discussed bv IWilson & Van Hammel (l2014h . 
who used AR Aur as an example of systems with a third com- 

pa nion. _ 

IWilson & Van Hammel (l2014t) also studied the W UMa type 
eclipsing binary XY Bootis as a prototype of a close double star 
with a nearly steady period change. The study clearly proves 
the advantage of simultaneous processing of all relevant data. 
That is why we want to compare the results of this sophisticated 
study describing period change of XY Boo during the time inter¬ 
val 1955-2009 with our results based on strictly the same obser¬ 
vational data. 


6.1. XY Bootis 

Eclipsing binary XY Bootis (= BDh- 20°2874 = HIP 67431; 
L max - 10.3 mag; Sp . E5V) was discovered as a variable star 
bv iHoffmeisteii (1 1 9351) . ITsesevichl (1 1 950i 1 1 9541) observed the star 
visually and classified it as an eclipsing binary of W UMa-type. 

The first p h otoelectric o bservations were obtained by 
iHindererl (Il960l) . I Wood! (1 19651) reanalysed Hinderer’s data and 
found the true period of the star P - tf37054. [Binnendii5 (Il97ll) 
observed six minima of XY Boo in BV, improved the ephemeris 
(Mo = 2440389.7321, P = 0?37055), and revealed remarkable 
increase in the period. IWinklej (1 19771) observed XY Boo also 
in 1976 in BV and derived thr e e othe r times of minima (see 
Table|6]l. lAwadalla & Yamasakil (1 19841) gave two new times of 
light minima and confirmed period increase. The history of the 
XY Boo investigati on from its discovery t o 1998 is described 
more e xtensiyely in iMolrk & WolfI (1 19981) . iMcLean & Hilditc^ 
(1 19831) measured spectroscopically the radial yelocity of both 
components and estimated the mass ratio q - 0.1 6 + 0.04. _ 

Th e first detailed period study was done by iMolik & Wolf! 
(Il998l) who used 43 moments of minimum spanning the inter- 
yal 1944-98 and calculated the quadratic ephemeris. The found 



- 0.2 0 0.2 0.4 0.6 0.8 

quadratic phase 


Fig. 10. BV curves of XYBoo (see the list in Table|5}. The phase is 
plotted according to our quadratic ephemeris (see Table|4j. The areas 
of points are proportional to their weights. Open squares are normal 
points (each of them represents the mean of about 50 measurements). 
Grey lines are fits by our phenomenological model (see Eq. (H and 
TableO. AS and AV display residuals of BV magnitudes and normal 
points from the LC phenomenological model. The scale of residuals is 
two times larger than the measure of BV light curves. 


record-breaking period increase P = L67(5)xl0 5.3 s per 

century was explained by mass transfer of 1.34 x 10“^ M© yr“^ 
from the secondary to the pri mary component. The results 
were confirmed and improyed by lYang et ^ (l2005l) who found 
P- 1.71 l(6)xl0“®= 5.4 s per century on the basis of standard 
0- C analysis of 54 minima times . 

IWilson & Van Hammel (1201 4t) combined in their calculations 
of the quadratic ephemeris of XY Boo phase information hid ¬ 
den in BV LCs obtained by BinnendiikI (Il97ll): Winkled (Il977l); 
lAwadalla & Yamasakil (1 19841) . RV c urves of McLean & Hilditchl 
(11983 ). all times of minima listed in lYang et al. (2005) and other 
33 minima timings collected in Table 1 of their article cover¬ 
ing the interval 2005-09. Using all of those data they found 
the mean period increase P - L6348(8)xl0“^ = 5.159(24) sper 
century. They also noticed some oscillation from JD 2 448 000 
to 2455 000 with no indication of periodicity. We have now 
collected many other observations proving the complexity of 
the XYBoo period variations. Nevertheless, for the sake of 
the comparison of the effectiveness of our method, we used 


in the following smal l study 
IWilson & Van Hammel (l2014l) . 


exclusively those data used by 


6.2. Phenomenological model of XYBoo variability 

We assumed, similarly as IWilson & Van Hammel(l2014l) . that the 
instantaneous period P(t) of XYBoo is lengthening with the 
constant rate P{t) - const. Then the phase function P{t) and the 
prediction of primary minimum times 0(E) can be approximated 
(according to Eqs. (|7]i, (I^) by simple relations: 

if 

= 0^,^Mo + PoE+\PaPE\ (17) 

2 ^ 
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Table 3. Nine parameters of phenomen ological (this paper) and syn¬ 
thetic LCs dWilson & Van Harnmdl2014h models. 


par. 

Model I 

W&VH 

prm 

solution I 

W&VH 

an [mag] 

0.102(4) 

0.098 

A 4 [mag] 

0.0064(6) 

0.0075 

a 12 [mag] 

0.112(4) 

0.104 

£>1 

0.0627(21) 

0.059 

a 2 l [mag] 

0.050(10) 

0.056 

D 2 

0.0558(20) 

0.053 

022 [mag] 

0.007(10) 

0.048 

r 

0.97(6) 

0.95 

A 3 [mag] 

0.1081(19) 

0.111 





For the meaning of parameters see Eq. jl 8 > . 

where t?i = (t-Mo)IPo, Mq is the JD timing of the basic primary 
minimum - the origin of counting of epochs E, Pq - P(t - Mq) 
is the instantaneous period at the time f = Mq. Integer doubling 
of epoch E is done (even for times of primary minima, odd for 
times of secondary minimafl 

Light curves of XYBoo with almost equally deep min¬ 
ima (see Fig.[IOll agree with W UMa classification. We 
found that all of them can be well fitted by the fol¬ 
lowing phenomenological model with nine free parameters: 
ail, an, 021 , 022 ,^ 31 ^ 4 , Di,D 2 , and P, 



phase 

Fig. 11. Models of B and V L Cs of XY Boo. Grey lin es are BV model 
light curves simulated by the IWilson & Van Hammd ( 120141) physical 
model, dotted lines correspond to our phenomenological model (see 
Eq. i fTSt and Table|4l(; dashed lines are hypothetical BV LCs of XY Boo 
without mutual eclipses. 

Table 4. Comparison o f common parameters derived by 
IWilson & Van Harnmd (l2014h and us. 


Flc{^) = ^{ak\ + akiA) ] 1 - ] 1 -exp 
k=\ V ^ 

-I- A 3 cos( 4 ; 7 r)?) -I- A 4 cos( 67 rj 5 -); (18) 

(pk^[P-(k- l)/2] - round [&-{k- l)/2], A = ^ - 1. 

deft 

The results achieved using this phenomenological model of 
XY Boo BV light curves are denoted as ‘Model I’. 

Other LC templates were BV synthetic LCs computed by the 
freely accessible W- D 2013 code with phy s ical p arameters of 
XY Boo published in IWilson & Van Hanimd (1201^ . The results 
obtained using these BV LC templates are denoted as ‘Model IF. 
Fig EI displays that the difference between the found phe¬ 
nomenological and synthetic BV model LCs are only marginal; 
numerically it is also documented in Table|3] 

For the description of shapes of radial velocity {RV) phase 
curves we used only a simple sinusoidal model (see Ea.(fT9T l) 
neglecting the Rossiter-McLaughlin effect (iMcLaughlinl 1 9241) 

RVki'ff) -Vy+As{k-\-q) sin(27r)?), (19) 

where k - \ refers to the first component, whilst k - 2x.o the 
secondary one, Vy is so-called y velocity, A 5 is the amplitude of 
the difference in the radial velocity between the components, all 
in km s“*, and q - OT 2 /'«i is the ratio of component masses. 

This approximation is fu lly sufficient for RV observations 
of iMcLean & Hilditchl (Il983h that cover only about 15 % of the 
phase curve. Results are in the bottom part of TableH) 

6.3. Used data and their weights. Modeis i and it 

Our period analyses of XYBoo was based on 1770 individual 
data points (see the list in Tabled acquired by three different 
techniques divided into 9 groups with various scatters cri and ctu 

® IWilson & Van Hammd (120141) used for modelling of phase function 
I? and its inversion 0 an exact solution of the basic equation, Eq. 0: 
&(l) = p-Un[l+P/Po(t-Mo)\, 0(E) = Po/P^p(EP- 1) - l] + Mo. 
Nevertheless, our models for & and 0 (see Eas. l l7b do not differ by more 
than 1.9 X 10“^ in the phase function and 0.3 s in time prediction from 
the exact ones. So they can be considered to be identical. 


par. 

W&VH 

Model I 

Model II 

all 

P 

1.6348(77) X 10^'’ 

1.6328(52) X 10^'' 

1.6442(47) X 10^'' 

Po 

O?370 560 205(18) 

Of370 560 213(14) 

O?370 560 198(12) 

M* 

0.35163(28) 

0.35185(26) 

0.35143(27) 

LC 

P 

1.497(37) X 10^'^ 

1.506(24) X 10-'^ 

1.502(25) X 10-'^ 

Po 

0f370559 889(92) 

0f370 559 913(70) 

0f370559 898(71) 

m; 

0.35162(29) 

0.35174(30) 

0.35164(30) 

Tm 

W&VH 

this paper 


P 

1.729(27) X 10^'^ 

1.745(26) X 10^'^ 


Pq 

0?370 559 58(11) 

0?370559 86(10) 


M* 

0.3565(21) 

0.3490(11) 


RV 

W&VH 

this paper 


Vy 

7.7 ± 6.6 

8±7 


4 

0.159 ±0.048 

0.154 ±0.037 


As 

286 ± 24 

282 ± 23 



Mq = HJD — 2 444 716, Vy and A 5 are expressed in km s * 


with respect to phenomenological (I), and synthetic (II), model 
phase curves. 

Scatter cti used to be the same or a bit smaller than the 
scatter crn with respect to c omputed synthetic phase curves 
(IWilson & Van Hammel 120141) . Scatters of radial velocity mea¬ 
surements - ctrvi and cr^^yn 28 and 30 km s“' - seems to be a 
good estimate of their values. The mean uncertainty of eclipse 
time^ cr-pnjin with respect to the model Eq. (fTTl i is (f?0043 = 
6.2 min. It should be mentioned that the inner uncertainty of in¬ 
dividual times of minima is typically one order better. The un¬ 
certainty cr-i-min is caused mainly by erratic fluctuations in the 
rate of the orbital period and the shape of LCs. The character of 
these changes is well displayed in Fig. [121 

Scatter ctj n for an individual subset of data (see Table|5]) 
served as the basis for weighting of each used measurement and 
interconnection of data of various nature. It enabled applying of 
regression computation outlined in Sect.|5] The computation 
represents the solution of 21 (case I) or 12 (case II) non-linear 
equations of 21 or 12 free parameters. The difference in the 
number of parameters needed is a result of the fact that in case 

^ We did not distinguish between primary and secondary light min¬ 
ima because their depths are nearly the same. 


1 - cosh 
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Fig. 12. XY Boo timing residuals from the quadratic ephemeris (with a 
P term). Published minima timings are denoted as squares, timings de¬ 
rived from LCs are circles which areas are proportional to their weight. 

Table 5. The list of used observational data including their source, 
number, specification, and scatter with respect to phenomenological 
model phase curves (Model 1) and synthetic phase curve (Model 11). 


source 

type 

N 

spec 

CTl 

O-II 

Binnendiik (19^) 

light 

440 

B 

0.011 mag 

0.011 mag 


curves 

441 

V 

0.0084 mag 

0.0087 mag 

Wmkler C1977_) 

light 

99 

B 

0.021 mag 

0.020 mag 


curves 

109 

V 

0.022 mag 

0.023 mag 

Awadala (1984) 

light 

280 

B 

0.021 mag 

0.023 mag 


curves 

284 

V 

0.025 mag 

0.026 mag 

McLean (1983) 

radial 

12 

I 

28kms-‘ 

28kms-' 


vel. 

10 

II 

30 kms“* 

30 km s“^ 

Yana et al. (2005) 

-hW&V(2014) 

eclipse 

timings 

75 

I+II 

0.0043 d 

0.0043 d 


Fig. 13. XY Boo timing residuals from the linear ephemeris (without a 
P term). Published minima timings are denoted as squares, timings de¬ 
rived from LCs are circles which areas are proportional to their weight, 
and a timing derived from radial velocity curves is noted as a diamond. 


Table 6. Comparison of the eclipse timings of XY Boo. 


publ. timing 

Model I 

Model II 

(0-C)i[d] 

epoch 

aut 

39 950.8121 

.81207(16) 

.81207(16) 

-0.00024 

-12860.5 

B 

39 951.9243 

.92427(14) 

.92430(13) 

0.00030 

-12857.5 

B 

39 953.7763 

.77635(13) 

.77628(13) 

-0.00038 

-12852.5 

B 

39 953.9626 

.96259(20) 

.96261(21) 

0.00058 

-12852.0 

B 

40 298.9470 

.94746(16) 

.94745(16) 

0.00089 

-11921.0 

B 

40 389.7319 

.73142(13) 

.73145(13) 

-0.00065 

-11676.0 

B 

42 569.7091 

.70875(55) 

.70875(55) 

0.00212 

-5793.0 

W 

42 577.6745 

.67530(65) 

.67535(65) 

0.00173 

-5771.5 

W 

42 582.6769 

.6750(12) 

.6750(12) 

-0.00108 

-5758.0 

w 

45 131.3793 

.3800(11) 

.3800(10) 

0.00047 

1120.0 

A 

45 132.3056 

.30470(80) 

.30450(85) 

-0.00126 

1122.5 

A 

44 325. 

.446(16) 

.446(16) 

0.03500 

-1055.0 

RV 


II we adopted template LCs from paper fwison & Van Hamrii^ 
(l2QI4l) . This means that we fixed all 9 parameters that we needed 
for the determination of BV template LCs. 

In the first part of solutions I and II we obtained the param¬ 
eters of the unified XYBoo ephemerides; Mq, Pq, and P (see 
Table|4]i, wh ich can be compared with the s ame parameters listed 
in Table 3 of lWilson & Van Hammel (1201 4h paper. Computed pa¬ 
rameters of the B, V, and PV phase curves according to these 
ephemerides are in Table[3]and Table|4] respectively. In the sec¬ 
ond part of the solutions we derived 12 virtual eclipse times and 
their uncertainties by means of the model Eq. (0 from BV pho¬ 
tometry and BV data - see Tablej^ Fig. [T^ [13] 


6.4. Discussion of resuits. Comparison with W&VH soiution 
6.4.1. General remarks 


We discuss here the results of three methods for period analy¬ 
sis that enables the processing of various type of data containing 
phase information. For their consistent comparison we used ex- 
actly the same observa t ional data set specified in the paper of 
IWilson & Van Hammel (l2014l) and the assumption that the pe¬ 
riod P rises uniformly with time: P(t) - Pq + P (t - Mq), where 
P - const. 

Model I uses the phenomenological model for the BV LCs 
(E g. (ITSll), whilst hybrid Model I I uses synthetic LCs computed 
bv lWilson & Van Hammel (1201 4l). This tec hnique is the one used 
e.g. bv izasche et aiT ( 2014h : ^^hel (l2015h . Both models assume 


Published timings of [ Binnendiild Il971h . IWinkleJ Il977h . and 
lAwadalla & Yamasakil lull are compared with our results (only the 
fraction is given) and their uncertainties derived from LCs and RV curves. 
(O-C)i is the residual of the found eclipse timing to quadratic ephemeris 
according to the Model II, N is the number of photoelectric observations used 
for the timing determination. 


strictly circular orbits of components and synchronous rotation. 
They also suppose that the shape of LCs given, namely, by the 
variable geometry of the orbiting system, are more or less con¬ 
stant. 

On the contrary, detailed inspection of the LCs proves that 
there are apparent changes on v arious time scales (see Fig.fTTI 
and BV light curves in Fig. 3 in I Wilson & Van Hammel 12014 ). 
Some of them can be attributed to instrumental effects like red 
noise (iPont et al.ll2006h or incomplete detrending of LCs. Others 
are due to intrinsic changes of the object itself, such a chromo¬ 
spheric and spot activity or shifts due to unsteady mass transfer 
between components (see the right part of Fig.lTSli. 


6.4.2. Unified ephemerides. Virtual eclipse times 

All quantitative results achieved by the three discussed models 
which can be directly compared are listed in Table|4] It is appar¬ 
ent that especially the increase of the orbital period P, the period 
P() - P{t = Mq), and the time of the basic primary minimum Mq 
derived from all data, from LCs only and from eclipse times, are 
identical within their uncertainties. The results derived from the 
analysis of radial velocity measurements are in very good con- 
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Fig. 1 4 . Residua of BV light curves of iBinnendiikI (Il97lh : Iwinkled 
l ll977i) : lAwadalla & Yamasakil ( 119841) clearly show variations in LC 
shapes. B and V residuals are denoted as circles and open squares, re¬ 
spectively. 


cordance as well. We are convinced that this conclusion is valid 
not only for XYBoo but also for other eclipsing systems. 

The re is an apparent identical disc repancy found by this pa¬ 
per and IWilson & Van Hammel (l2014l) in the value of the mean 
deceleration P derived only from LC shifts and/or eclipse times 
(Tm) (see Table|4] rows LC and T^n): Flc = 1.506(24) x 10“^, 
= 1.745(26) X 10^®. 

The first possible explanation, that it is the result of different 
distribution of LCs and times data along the 0-C diagram (see 
Fig . [T3l [12]) and the inconstancy of P where P > 0), has proven 
to be false. We tested the period models with P term and found 
that P - 0(4) X 10“'^ d“'. The effect is then v ery likely caused 
by the presence of the pair of eclipse minima of lHinde^ (1 19601) 
at the very beginning of the 0-C diagram (see Fig. (Ell that sig¬ 
nificantly deviates from its parabolic prediction. However, our 
thorough inspection of the original Hinderer’s measurements 
showed that both eclipse times are correct. 

Models I and II enabled calculation of ‘virtual eclipse times’ 
for selected subsets of photometric and radial velocity data us¬ 
ing the relation Eq. (|9]). These times, including their uncertain¬ 
ties, are listed in Table|6| together with eclipse t imefl p ublished 
by authors of BV photometries -iBinnendii^ (119711) : IWinkled 
(ll977h : [Awadalla & Yamasal3(ll984l) The results as well as their 
uncertainties obtained by Model I and Model II are nearly iden¬ 
tical. The last ones are very important as they enable the evalu¬ 
ation of real inner accurac y of individual times, th e construction 
of precise 0-C diagrams (iMikulasek et alJl201 lah . and the dis¬ 
cussion of subtle problems of stability of the period, variability 
of LCs etc. There is no systematic difference between the pub¬ 
lished eclipse timings (see the first column of Table|6ll and those 
derived by our Models. 

The excellent agreement of comparable results obtained by 
all three models proves that they can be used as good alternatives 
in solving common tasks of the eclipse system period analysis. 


All the eclipse times were given without quoting their uncertainties. 


6.4.3. BV an(d RV phase curves 

The differences between the phenomenological and the synthetic 
BV light curves are insignificant - they represent about 0.5% of 
the amplitude of the light changes (see Fig.ITTI). The parame¬ 
ters of the phenomenological model of both mentioned LCs are 
in Table|3 It seems that phenomenological modelling is a good 
method for expressing common types of LCs using only a very 
small number of parameters. 

The fit of the observed BV light curves by Model I is fairly 
good (see Fig. doji - scatter of normal points in B and V are bet¬ 
ter than 0.004 and 0.003 mag. The larger part of this scatter is 
caused by the erratic variability of LCs due to stellar activity 
and inconstant O’Connell effect. The changes of LC shapes are 
clearly visible in the phase diagram of BV residuals (Fig. [T4l) for 
all three photometric data sets. _ 

Th e detailed inspection of measurements of iBinnendii^ 
(Il97lh shows weak and variable O’Connell effect (data from 
1968 and 1969 disp layed O’Connel l effe ct of the opposite sign), 
whilst residuals of IWinklej (Il977h and lAwadalla & Yamasakil 
(1 19841) are dominated by double waye going in antiphase. It 
seems that such seasonal yariations in light curye shapes are 
quite common. Sometimes they can be much more dramatic (se e 
e.g. changes in LCs of notorious XY UMa. iLister et al.ll200lT) . 
These LC changes lowered the accuracy of the fit and eclipse 
timing determination. 

Yet unp ublished obseryations of XY Boo obtained during 
one month (IMikulasek et al.ll2016l) proye that the time scale of 
LC changes is usually shorter than seyeral days. 


7. Conclusions 

Phenomenological modelling presents an admissible alternatiye 
for the solution of selected ES research tasks. We state that: 


1 . 


2 . 


3. 


4. 


5. 


The estimation of parameters and their uncertainties ob¬ 
tained by our phenomenological (and also hybrid) modelling 
and other well-proyen methods, including solutions by so¬ 
phisticated physical models of eclipsing systems, are almost 
the same (see e.g. Sect. l4.3l and l6.41 i. 

Phenomenological modelling is based on the minimization 
of the sum, which enables the simultaneous processing 
of different sources of phase information (complete LCs and 
their segments, radial yelocity curyes, indiyidual mid-eclipse 
times etc.). 

Simple model function of eclipses (Eq. ([T3l l) can be also used 
for good determination of mid-eclipse times and their un¬ 
certainties of indiyidual obseryations of stellar eclipses (our 
model has been standardly used for the determination of 
time minima of original obseryations of EBs by Variable star 
and E xoplanet Section of Czech Astronomical Society since 
2011, iBrat et al.ll2012l) 


As a by-product of the phenomenological modelling we ob¬ 
tain the list of yirtual minima times deriyed from LC data 
which help us i.a. to quickly check the selected phenomeno¬ 
logical model of LCs or period yariability (see Sect. l6.4.2l l. 
Light curye model parameters can be used for the apt de¬ 
scription of both one-colour and multicolour LCs of eclips¬ 
ing systems. Eor example, shapes of ten curyes of HST 
spectrophotometry (320-9 80 nm) of the exoplan et transit of 
HD 209458b, taken from iKnutson et a k ( 20071). are deter¬ 
mined by only nine parameters (see IMikulasek et alJl201^ . 
We can use this application e.g. for description and classifi¬ 
cation of obseryed LCs of eclipsing systems. 
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6. Knowing the template of LCs (from physical models or ob¬ 
served LC of superior quality) we can quickly modify the 
presented phenomenological model and establish a hybrid 
model (Model II in Sect.|6]l with the diminished number 
of free parameters. As we showed in Sect.|6] the results of 
this approach are the same results as physical or pure phe¬ 
nomenological modelling. However, application of the latter 
method is much faster. 

7. Phenomenological (and also hybrid) modelling of an ES 
could solve standard tasks of ES research (based on every 
sources of phase information), especially the improvement 
of ES ephemeris for standard period models of; 

o systems with constant period, 

o systems in a steady regime of mass and angular momen¬ 
tum transfer between components (e.g. XY Boo, Sect.|^, 
o systems with erratic changes of period - phase function 
is here approximated using so-called 0-C time shifts of 
the observed phase curves versus the predicted LC de¬ 
rived by the period model with fixed parameters (typi¬ 
cally Pq, Mq or Pq, Mq, P) - see Eq. (|9]), and Sect. l6.4.2l 
o eclipsing systems influenced by the gravitational attrac¬ 
tion of a third body (see the pre liminary paper about 
AR Aur in iMikulasek et al.11^1 Ibi and the discussion in 
Wilson & VanHamme 2014jJ 


Phenomenological modelling of eclipsing systems presented 
above has also its disadvantages and limitations: 


The method does not provide direct information about the 
geometry and physics of eclipsing systems, especially the 
inclination of the system, relative radii, temperatures and the 
form of components. 

The phenomenological models of EBs light curves in their 
proposed form are not able to achieve an accuracy better 
than 0.5 % of their total amplitude of them (the accuracy of 
Kepler or CoRoT light curves is better by one or more or¬ 
ders). Nevertheless, such uncertainty is fully acceptable for 
many application. 

It may be a bit troublesome for beginners to select the opti¬ 
mum model of LCs. Using an inappropriate model may lead 
to nonconvergent or bad solutions. 

Phenomenological modelling of ESs has not been made 
available in such a user-friendly format as other modern 
physical models (e. g. Hadra^l2004iPrsa et al]|2008ll20Tlt 


Iwilson & Devinr^IlQTlh until 


Phenomenological modelling, influenced by modern phys- 
ical models, has be en developed gradually since 2008 
jMikulasek et al.ll2008h . its elements and principles were used 
in the majority of papers of this paper’s author. The detailed ap¬ 
plications of the presented bases of the method will be published 
elsewhere. 
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